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Abstract 
This study explores the uncertainties associated with general circulation 


model (GCM), emission scenarios and natural climate variability into 
estimates of climate change impacts over two Mexican catchments. The 
two selected catchments have contrasted climate patterns and 
hydrological regimes. Climate ensemble simulations from three GCMs 
under three SRES scenarios (B1, A1B, and A2) are considered to feed 
the physically-based semi-distributed SWAT hydrological model. The 
simulations cover the 30-year reference period (1971-2000) and two 
20-year periods (2046-2065 and 2081-2100) in the future. For the set 
of hydrological indicators, the results show that the high flows are more 
prone to be influenced by the GCM uncertainty. A weak sensitivity of the 
hydrological indicators is observed to emission scenarios over the study 
catchments. We find that the uncertainty related to the natural climate 
variability should not be neglected in climate change impact studies, and 
appears, to some extent, to be as critical as the GCM and scenario 


uncertainties for high flows. 


Keywords: Climate change impact, model chain, streamflow, 


uncertainty, natural variability, SWAT, Mexico. 


Resumen 
Este estudio investiga la incertidumbre asociada al modelo de circulación 


general (MCG), al escenario de emisión y a la variabilidad natural del 
clima en la estimación del impacto del cambio climático en dos cuencas 
mexicanas. Las cuencas seleccionadas tienen diferentes climas y 
regímenes hidrológicos. Las simulaciones climáticas de tres MCG bajo 


tres escenarios de emisión (B1, A1B y A2) son usados para alimentar el 
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modelo hidrológico SWAT el cual es un modelo del tipo físico y 
semidistribuido. Las simulaciones cubren 30 años en un periodo de 
referencia (1971-2000) y dos periodos de 20 años en el futuro (2046- 
2065 y 2081-2100). Los resultados muestran que los caudales altos son 
más propensos a ser influenciados por la incertidumbre asociada al 
MCG, además, los indicadores hidrológicos tienen una sensibilidad 
menor al escenario de emisión para las cuencas de estudio. Finalmente, 
la incertidumbre asociada a la variabilidad natural no debe pasarse por 
alto en los estudios de cambio climático ya que, hasta cierto punto y 
para caudales altos, esta puede ser tan crítica como la incertidumbre 


asociada al MCG y al escenario de emisión. 


Palabras clave: impacto del cambio climático; cadena de modelación; 


escurrimiento; incertidumbre; variabilidad natural; SWAT; México. 
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Introduction 


Climate change projections for Mexico and Central America suggest 
changes in precipitation and an increase in temperature that will affect 


the future availability of water resources in the region (IPCC, 2014). 
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Table 1 shows some studies of climate change impacts on water 
resources in Mexican basins. The use of several General Climate Models 
(GCMs) allows the uncertainty related to the climate model in the 
evaluation of climate change impacts on water resources to be 
assessed. Recently Velázquez-Zapata, Troin and Dávila-Ortiz (2017), 
and Molina-Navarro et al. (2016) evaluated climate change impacts with 
five and six climate models, respectively. On the other hand, the work of 
Rivas, Gúitrón and Montero (2011) combined outputs from 23 GCMs to 
obtain an averaged climate model simulation. However, Arnell and 
Gosling (2013) claimed that averaging model outputs lead to a loss of 
signal in precipitation (when estimated changes are of different sign) 
and to an underestimation of the spatial variability of the signal. 
Furthermore, none of the previous studies have considered the 
uncertainty associated with the natural climate variability (NCV). The 
NCV is estimated by repeating a climate change experiment using a 
given GCM several times, when only the initial conditions are changed 
by small perturbations (Braun, Caya, Frigon, € Slivitzky, 2012). In that 
aspect, Deser, Knutti, Solomon and Phillips (2012) showed that the 
natural fluctuation of climate contributes substantial uncertainty in the 
projections of climate change over North America, especially on high 
latitudes. In contrast, the authors showed that Mexican locations are 
less subject to the NCV. 


Table 1. Studies of climate change impact on water resources over Mexican basins. 
The table shows the number of Global Climate Models (GCMs) and Regional Climate 
Models (RCMs) used in the study. 


Reference Basins Indicators Climate models Emission Horizon Natural 
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Scenario Variab. 
Mendoza The whole Aridity index 2 GCMs No No 
Villanueva country vulnerability 2050 
and Adem divided in 12 (e.g., water 
(1997) hydrological availability 
regions and demand) 
Annual 
surface runoff 
Gratiot et al. | Cointzio (650 Aridity index 1 GCMs A2 Decades No 
(2010) km?) centered 
in 2030, 
2060 
and 2090 
Rivas, Lerma- annual runoff Average of 23 A1B, A2 2030 and No 
Guitrón and Chapala (54 volume GCMs 2050 
Montero 450 km?) 
(2011) 
Maderey, Lerma- 2 GCMs and a The scenario 2025- No 
Jiménez and Chapala- Avalaibla thermodynamic considers 2050 
Carrillo Santiago Suma Wáter model duplication of 
(2013) (135 836 CO, content in 
km2); Balsas the 
(117 638 atmosphere. 
km?), and 
Pánuco (98 
302 km?) 
Tapia, Yaqui River Water 1 GCM A1B and A2 2010- No 
Minjarez, Basin (74 balance 2099 
Espinoza 054 km?) 
and Minjarez 
(2014) 
Velázquez, Tampaón Mean monthly 1 CGM-RCM A2 2071- No 
Troin and River Basin streamflow 2100 
Caya (2015) | (23 373 km?) indicators 
Robles- Sonora River Reservoir 1 GCM-RCM A2 2031- No 
Morua Che, Basin (20 inflows 2040 
Mayer and 648 km?) 
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Vivoni 
(2015) 
Molina- Guadalupe Water 6 GCMs Bi and A2 2010- No 
Navarro et River Basin balance 2039 and 
al. (2016) (2380 km?) 2070- 
2099 
Velázquez- Tampaón Mean monthly 5 GCMs A2 2071- No 
Zapata et al. River Basin streamflow 2100 
(2017) (23 373 km?) indicators 


The uncertainty associated with the NCV had been evaluated in 
several studies of climate change impacts on water resources in higher 
latitudes; (e.g., Velázquez et al., 2013). However, to our knowledge, the 
uncertainty associated to the NCV on hydrological impacts of climate 


change has not been evaluated over Mexican basins (Table 1). 


The aim of this study is to assess the uncertainties derived from 
the climate change impacts on the hydrological regime of two contrasted 
Mexican catchments. Simulations from three GCMs under three SRES 
(Special Report on Emissions Scenarios) (IPCC, 2010) scenarios (A2, Bl, 
and A1B) are used to feed the physical-based semi-distributed SWAT 
(Soil and Water Assessment Tool) hydrological model over a reference 
period (1971-2000) and two future periods (2046-2065 and 2081- 
2100). Three sources of uncertainty are compared in this study: the 


GCMs, the SRES emission scenarios and the NCV uncertainties. 


Of Table 1, we can see that most impact studies evaluated 
changes on the basins” water balance. In the present study, the basins” 
hydrology is analyzed through two specific indicators: annual mean 


flows and high flows. 
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The manuscript is organized as follows: section 2 presents the 
study area, the observational and climate model datasets; section 3 
gives a description of the hydrological model and the bias correction 
procedure; section 4 presents significant results on the evaluation of 
climate change impacts on the basins” hydrology and the associated 


uncertainties; and section 5 concludes. 


Experimental design 
Study catchments 


The Tampaon River Basin is located in East-Central Mexico (23 373 km; 
IMTA, 2018), lying mostly in the State of San Luis Potosí. The Tampaon 
River Basin is a sub-catchment of the Panuco River Basin, which flows 
into the Gulf of Mexico. Its topographic relief ranges from 18 m to 3500 
m, with maximal elevation on the mountainous eastern watershed 
(Figure 1). The Tampaon River Basin has a variety of climatic regions, 
as a result of topographic variation in the Sierra Madre Oriental 
mountain chain (Sedue, 1986;): the farthest headwater of the Tampaon 
River is located in an arid cold steppe (BSk) climate type region and 
flows west to east through a temperate dry winter hot summer (Cwa) 


climate type region (Peel, Finlayson, 8£ McMahon, 2007). Then the Valles 
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River, which is located in tropical Savannah (Aw) and tropical monsoon 
(Am) climate type regions, joins the Tampaon River providing the 


largest volume of streamflow. 


Tampaon River Basin 


United States of America 


Elevation 
Max: 3 500 


Min : 18 
/5x, Meteorological stations 


Mexico (+) Gauging station 


Elevation : 
- Max: 3529 WR ad E 4 Po 
Meteorological stations 
_— - 5 Pacific Ocean 


(+) Gauging station 


Figure 1. Location of the study catchments. 


The Papagayo River Basin is entirely located in the state of 
Guerrero (South Mexico), covering 7 067 km2 (IMTA, 2018), with a 
maximum elevation of 3529 m (Figure 1). The Papagayo River Basin is 


located in a Tropical Savannah (Aw; Peel et al., 2007) climate type 
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region, and it flows from the Sierra Madre del Sur headwaters to the 


Pacific Ocean. 


Observational dataset 


The daily time series of minimum and maximum temperatures and 
precipitation were taken from the CLICOM climatologic database, built 
by the SMN (Servicio Meteorológico Nacional) via the CICESE (2017) 
website. The data come from seven stations for the Tampaon River 
Basin and from five stations for the Papagayo River Basin over the 
1971-2000 period (Figure 1). 

The mean annual temperature is about 21.5% C and 24.1% C for 
the Tampaon River and the Papagayo River, respectively. The coldest 
month is January and the warmest month is May (Figure 2a). The rainy 
season extends from May to October (Figure 2b) for both catchments. 
The area's mean annual rainfall is about 1080 mm and 1540 mm for the 


Tampaon River and the Papagayo River, respectively. 
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a) mean monthly temperature b) mean monthly precipitation 
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Papagayo 


La Ll mM 
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month month 


Figure 2. Mean monthly temperature and precipitation over the study 


basins. 


The discharge data were obtained from the National Database of 
Surface Water (í/.e., BANDAS dataset; IMTA, 2018) for the 1971-2000 
period. The observed mean monthly discharges at the gauging stations 
over the two basins is showed in Figure 3. The Tampaon River Basin 
presents two peak flows in July and September, while the Papagayo 
River has its maximum flow in September. For both catchments, the 
September's peak flow is about 450 m* s* (47 mm month? and 164 
mm month'* in normalized-area streamflow for the Tampaon River Basin 


and the Papagayo River Basin, respectively). 
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a) mean streamflow b) mean streamflow 
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Figure 3. Observed mean monthly discharges for the study catchments over the 
1971-2000 period. 


The climate model ensemble simulations 


The climate model ensemble simulations were derived from three GCMs 
through the CMIP3 multi-model dataset under the SRES three emission 
scenarios (A1B, A2 and B1; Meehl et al., 2007; Table 2). The GCMs and 
SRES emission scenarios are considered to be equally plausible in order 
to evaluate the uncertainty of the impacts of climate change on water 


resources over the study catchments. The GCMs were selected because 
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they provide several simulations for each emission scenario, which are 


used to assess the uncertainty related to the NCV., 


Table 2. CMIP3 models used in the study (see Meehl et al. (2007) for full references). 


Institute CMIP3 1.D. study 1.D. members 


Meteorological Institute of the University of Bonn, 


Meteorological Research Institute of KMA, 


and Model and Data group ECHO-G ECH 3 

Canadian Centre for Climate Modelling € Analysis CGCM3.1 CcGC 5 
MRI- 

Meteorological Research Institute CGCM2.3.2 MRI 5 


The use of direct GCM outputs as input to hydrological models can 
lead to unrealistic streamflow simulations due to the systematic errors 
(biases) of climate models (Jung, 2005). Therefore, before the use of 
GCM's outputs in hydrological models, it is necessary to adjust climate 
models' distribution in order to fit the observations” distributions through 
statistical post-processing methods (Teutschbein, Wetterhall, 8 Seibert, 
2011). 


In the present study, the method used for the post-processing and 
downscaling of the climate model outputs is the daily translation (DT) 
method (Mpelasoka 8: Chiew, 2009). In DT method, the transfer function 
is obtained with the differences (in percentiles) between the frequency 
distributions of observations and GCM outputs in the historical period. 
Then, those differences are applied to the GCM frequency distribution in 


the future period. When DT method is applied on a model with an 
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ensemble of simulations, the transfer function is computed on the 
ensemble mean and then each member is adjusted separately. The 


corrected meteorological variables are corrected in reference period as: 


ref _ rpref ref ref 
Iori) a Tom d) + a ma) Em ed (1) 


F F pref 
re _ pre obs(m,q) 
Eorrid) — ¿sim(d) (E ) (2) 


where Tcorr Y Pcorr are the corrected temperature and precipitation in 
reference (ref) period for a given day (d), and m, q, obs and sim are 
month, percentile, observations and climate simulations respectively. In 


future (fut) period, the meteorological variables are corrected with: 


fut _ rpfut ref ref 
Teorr(d) — ¿sim(d) + La ma) as a (3) 
pfut  _ pfut KR od 4 
corr(d) "— *sim(d) a ( ) 
sim(m,q 
Methods 
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The hydroclimate model chain is based on an ensemble of bias-corrected 
climate simulations from three GCMs under the three SRES scenarios 
(B1, A1B and A2). Those simulations are used to feed the hydrological 
model. The resulted hydrological simulations are analyzed to evaluate 
the uncertainties related with the GCM, the emission scenario and the 
NCV into the estimates of climate change impacts on the hydrological 


indicators. 


The hydrological model 


The Soil Water Assessment Tool (SWAT) (Arnold, Srinivasan, Muttiah, € 
Williams, 1998) hydrological catcthment model was selected in this 
study. SWAT is a semi-distributed physically-based model extensively 
used to evaluate the climate change impacts on water resources (e.g., 
Githui, Gitau, Mutua, Bauwens, 2009; Troin, Velázquez, Caya, 8 
Brissette, 2015). SWAT distributes the catchment in several Hydrologic 
Response Units (HRUs) in order to take into account the spatial 
variability (topography, land use, and soil type) of the catchment. For 
each HRU, SWAT simulates the basin's water balance in a daily step 
based on daily meteorological variables. The water balance is expressed 
as: 
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W, =W, + Diz (Pi — Qisup — ET¡ — W¡ — Qigw) (5) 


where Wt is the soil moisture content at time t; W, is the initial 
soil moisture content. For a given day /, P, is the precipitation, Qisurr iS 
the surface runoff, ET; is the evapotranspiration, w; is the water 
percolated through the soil profile and Qigw is the groundwater flow. All 
terms are expressed in mm of water. Further details of SWAT's water 
balance computation are given in Neitsch, Arnold, Kiniry, Williams y King 
(2002). 


For both basins, SWAT was calibrated over the 1971-1985 period and 
validated over the 1986-2000 period. The performance of the 
hydrological model was assessed with the Nash-Sutcliffe (NS) 
coefficient. Overall, SWAT performs quite well at simulating streamflow 
over the calibration and validation periods, with NS values of 0.91 and 
0.85 for the Tampaon River Basin, respectively, and of 0.86 and 0.84 for 


the Papagayo River Basin, respectively. 


Hydrological indicators 


Two hydrological indicators were selected to evaluate the climate 


change impacts on the catchments” hydrology: 
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1. The mean daily runoff (overall mean flow; OMF) 


2. The flow exceeded on average every two years (HF2). To compute the 
high flow indicator it is supposed that the streamflow series follow the 
log Pearson III PDF (e.g., Muerth et a/., 2013; Troin et al., 2015). 


The GCMs used in this study provide an ensemble of climate 
simulations which are used to evaluate the uncertainty related to the 
NCV. In order to achieve this task, the relative changes on the 
hydrological indicators are computed with permutations, as in Velázquez 
et al. (2013). This method assumes that each member of the ensemble 
(in both the reference and future periods) is an independent proxy of 
climate. Therefore, the permutation allows to compare the future of a 


given member with the present of all members. 


The relative differences of the hydrological indicators (Aj) are 


computed as follow: 


a 


A 


ij per (6) 
] 


In this equation, the / and j subscripts represent the climate model 
simulations used to compute the hydrological indicators. Hence, CGC 
and MRI obtain 25 values of relative differences while ECH obtains 9 
values. The boxplot of relative changes (as in Figure 4) gives an 


estimation of the uncertainty related to the NCV., 
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NS 


Change 


ij permutations 


Figure 4. Permutation process” diagram. 


Climate change's impact on  hydrological indicators (AI) is 
expressed as the difference between simulated indicators from the 


reference (1"2f) to the future period (14%), 


Results 
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Future changes in precipitation and temperature 


The changes in the climatological variables for the future period are 
estimated by analyzing the differences in the bias-corrected GCM- 
simulated precipitation and temperature between the future (2046-2065 
and 2081-2100) and reference (1971-2000) periods. Figure 5 and 
Figure 6 show the climate change signals for mean annual temperature 
and precipitation over the Tampaon River and the Papagayo River 
Basins, respectively. The results indicate that, for both future periods, 
the temperature is expected to increase; however, the GCMs predict 
different values of increased temperature for a given emission scenario. 
The three emission scenarios used in this study (¡.e., A1B, B1 and A2) 
consider different assumptions of future global population, gross world 
product and energy production (IPCC, 2010). As expected, the largest 
changes in temperature are predicted for the horizon 2081-2100 under 
the A2 emission scenario, following by the A1B and B1 scenarios. The 
increase in temperature also varies depending of the GCM: ECH and 
CGC generally present the smallest and largest increase in temperature, 


respectively. 
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b) A1B (Horizon 2046-65) 
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Figure 5. Climate change signals over the Tampaon River Basin as estimated by the 
three GCMs for the horizons 2046-2065 (upper panels) and 2081-2100 (lower panels) 
under the Bi(a and d), A1B (b and e) and A2 (c and f) scenarios. 
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Figure 6. Same as Figure 5 but for the Papagayo River Basin. 
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Results are somewhat more mixed for precipitation; from Figure 5 
and Figure 6 it could be seen that, over the two basins, CGC predicts a 
decrease in precipitation, while ECH predicts an increase in precipitation. 
On the other hand, MRI shows small changes in precipitation over the 


two basins. 


In addition, Figure 5 and Figure 6 show that the spread of the 
relative changes both in temperature and precipitation for the horizon 
2081-2100 is larger than the spread of relative changes for the horizon 
2046-65, especially for the A2 scenario. This highlights how climate 
model uncertainty varies when considering the prediction horizon and 
the emission scenario. Furthermore, it can be seen that the relative 
changes in precipitation and temperature obtained with the ensemble 
simulations can overlap the relative changes in precipitation and 
temperature of another GCM (e.g., horizon 2046-2065 under the B1 
scenario; Figure 5a), so that the NCV uncertainty is comparable to the 
GCM uncertainty (in this horizon under the B1 scenario). In contrast, the 
uncertainty related to the NCV becomes less important in the horizon 
2081-2100 (e.g., Figure 5f). This result is in agreement with the work of 
Hawkins and Sutton (2009), which analyzes different sources of 
uncertainty in the prediction of the surface air temperature. The authors 
showed that the relative weight of the main sources of uncertainty in 
climate impact studies (i.e., emission scenario, climate model and 
natural climate variability) varies with the prediction lead time. Hence, 
the dominant sources of uncertainty are associated with climate models 


and emission scenarios for prediction horizons of many decades. 
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Uncertainty associated with natural climate variability 
on hydrological indicators 


Figure 7 and Figure 8 show the relative changes on the selected 
hydrological indicators between the reference (1971-2000) and future 
(2046-2065) periods (as computed with Equation 6) over the Tampaon 


River Basin and the Papagayo River Basin, respectively. 
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Figure 7. Relative changes (as calculated with Eq. 6) of overall mean flow (OMF, 
upper panels) and 2-yr return period high flow (HF2, lower panels) between the 
reference (1971-00) and the future (horizon 2046-65) periods over the Tampaon River 


Basin. The boxplots depict the median and the 25th and 75th percentiles. 
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Figure 8. Same as Figure 7 but for the Papagayo River Basin. 


From these figures it could be seen that CGC generally leads to 
negative relative changes on OMF; however, ECH and MRI often provide 
positive and negative relative change values for this indicator. For 
instance, over the Papagayo River Basin under the A1B scenario (Figure 
8b), the simulated OMF values by ECH varies between -4.5% and 
+6.3%, with a median value of -0.1%. It can be noticed that, in some 
cases, the range of the relative differences of OMF obtained with one 
GCM covers the range of the relative differences obtained with another 
GCM. For example, over the Tampaon River Basin (A1B scenario; Figure 
7b), the relative difference of OMF simulated by MRI ranges from -5.5% 
to 17%, and covers the relative difference of OMF obtained with ECH 
(between 7% and 9%); however, the median relative differences of OMF 


are very close between both GCMs, with OMF values of about 8%. This 
22 


Tecnología y ciencias del agua, ISSN 2007-2422, 11(1), 01-36. DOI: 10.24850/j-tyca-2020-01-01 


2020, Instituto Mexicano de Tecnología del Agua 
Open Access bajo la licencia CC BY-NC-SA 4.0 
Tecnología y (https://creativecommons.org/licenses/by-nc-sa/4.0/) 


CienciaszAgua 


situation is more common for the high flow indicator (Figures 7d, 7e and 
7f), indicating that the climate change signals on HF2 obtained with 
different GCMs may be similar due to the NCV uncertainty. 


The Wilcoxon rank-sum test (Wilcoxon, 1945) is performed to 
assess if two samples obtained from two distinct GCMs are independent 
or not. The Wilcoxon test results showing the rejection or not of the null 
hypothesis (i.e., the two series are independent samples from identical 
continuous distributions with equal medians) are depicted in the 
following tables. Table 3 and Table 4 compares (in pairs), for a given 
GCM, the series of the climate change impacts on the indicators from 
the three emission scenarios (B1, A1B and A2) for the Tampaon River 
Basin and the Papagayo River Basin, respectively. Similarly, Table 5 and 
Table 6 compares (in pairs), for a given emission scenario, the series of 
the climate change impacts on the indicators from the three GCMs 
(CGC, ECH and MRI) for the Tampaon River Basin and the Papagayo 
River Basin, respectively. The gray areas in Tables 3 to 6 indicate a no 
rejection of the null hypothesis, that is, a similarity in the climate 


change signal due to the NCV. 


When emission scenarios are compared, the null hypothesis is not 
rejected in all cases for HF2 (Tables 3 and 4). Results for OMF are more 
mixed, since over the Tampaon River Basin for the 2081-2100 period, 
the null hypothesis is rejected for all cases for CGC but not for ECH and 
MRI (Table 3). 


Table 3. Wilcoxon test results (p-values) of the comparison of scenario pairs for the 
Tampaon River Basin. The gray areas indicate no rejection of the null hypothesis at 


significance level of 5%. 
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CGC ECH MRI 
Ind. Hor. B1- B1-A2 |A1B-A2 B1- B1-A2 |A1B-A2 B1- B1-A2 |A1B-A2 
A1B A1B A1B 


OMF |2046- 0.712 | 0.030 | 0.088 | 0.001 | 0.605 |<0.001 | <0.001 |<0.001| 0.229 
2065 


2081- |<0.001|<0.001| 0.002 |[<0.001|<0.001| 0.114 | 0.063 | 0.003 | 0.081 
2100 


HF2 [|2046- 0.954 | 0.222 | 0.174 | 0.796 | 0.436 0.863 | 0.125 | 0.222 | 0.907 
2065 


2081- 0.383 | 0.112 | 0.269 | 0.436 | 0.436 | 0.436 | 0.712 | 0.628 | 0.756 


2100 
Table 4. Same as table 3 but for the Papagayo River Basin. 
CGC ECH MRI 
Ind. Hor. B1- B1-A2 |A1B-A2 B1- B1-A2 ¡A1B-A2 B1- B1-A2 |A1B-A2 
A1B A1B A1B 


OMF |2046- |0.4263 | 0.0436 |<0.001| 0.077 [0.1359 <0.001 | 0.4263 | 0.0914 | <0.001 
2065 


2081- 0.415 | 0.017 |0.2443|<0.001| 0.014 | 0.6048 | 0.4492 | 0.0298 | 0.229 
2100 


HF2 [|2046- 
2065 0.6554 | 0.8461 | 0.6414 | 0.3401 0.0625 | 0.004 | 0.8766 0.4492 | 0.4263 


2081- 
2100 0.9227 | 0.393 | 0.3825 | 0.3401 0.3865 | 0.9314 | 0.6695 | 0.6004 | >0.99 


24 
Tecnología y ciencias del agua, ISSN 2007-2422, 11(1), 01-36. DOI: 10.24850/j-tyca-2020-01-01 


2020, Instituto Mexicano de Tecnología del Agua 
Open Access bajo la licencia CC BY-NC-SA 4.0 


Tecnología y > (https://creativecommons.org/licenses/by-nc-sa/4.0/) 
CienciaszAgua 


Table 5. Wilcoxon test results (p-values) of the comparison of GCM pairs for the 
Tampaon River Basin. The gray areas indicate no rejection of the null hypothesis at 
significance level of 5%. 


CGC-ECH CGC-MRI ECH-MRI 
Ind. Hor. B1 A1B A2 B1 A1B A2 B1 A1B A2 
2046- | <0.00 
2065 1 <0.001 | <0.001 | <0.001 | <0.001 | <0.001| 0.258 | 0.310 | 0.118 
OMF 
2081-| <0.00 
2100 1 <0.001 | <0.001 | <0.001 | <0.001 | <0.001 | <0.001 | 0.532 | 0.009 
2046- 
2065 | 0.086 | 0.029 | 0.019 | 0.214 | 0.015 | 0.001 | 0.292 | 0.458 | 0.435 
HF2 
2081- 
2100 | 0.001 | 0.001 |<0.001 | <0.001|<0.001|<0.001| 0.226 | 0.086 | 0.017 
Table 6. Same as Table 5 but for Papagayo River Basin. 
CGC-ECH CGC-MRI ECH-MRI 
Ind. Hor. B1 A1B A2 B1 A1B A2 B1 A1B A2 
2046- 
2065 <0.001 | <0.001 | <0.001 | 0.005 |<0.001 | 0.756 | 0.969 |<0.001 | 0.002 
OMF 


2081- 
2100 <0.001 <0.001 |<0.001 | 0.002 | 0.074 | 0.006 | <0.001 | <0.001 | <0.001 


2046- 
HF2 [2065 0.369 | 0.274 | 0.042 | 0.304 | 0.027 | 0.641 — 0.585 | 0.482 | 0.160 


2081- | 0.032 | 0.009 | 0.067 | 0.187 | 0.214 | 0.587 | 0.585 | 0.056 | 0.051 
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2100 


In addition, when GCMs are compared (Tables 5 and 6), results 
show a general rejection of the null hypothesis for OMF. The exception is 
made for the Tampaon River Basin at the 2046-2065 horizon, where the 
null hypothesis is not rejected for the ECH-MRI model pair in any 
scenario. For OMF over the Papagayo River Basin, the comparison of 
OMF series leads to no rejection of the null hypothesis in the three 
isolated cases. However, the results are different for HF2: over the 
Tampaon River Basin (Table 5), the null hypothesis is not rejected in 
most cases, when the ECH-MRI model pair is tested. Similarly for the 
Papagayo River Basin (Table 6), most of the CGC-MRI model pair and all 
the ECH-MRI model pair lead to a no rejection of the null hypothesis. 


These results show a large uncertainty on the high flow indicator, 
but this uncertainty is more related to the NCV than to the choice of the 
GCM or the emission scenario. On the other hand, we can see that the 


NCV uncertainty is less important for OMF than for HF2. 


Conclusion 


The scope of this work ¡is the evaluation of the impacts of climate change 


on the hydrology of two contrasted catchments in Mexico, by 
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considering three types of uncertainty: climate model, emission scenario 
and Natural Climate variability (NCV). 


Regarding the projected changes on meteorological variables, the 
analysis shows a large variability on precipitation over the two 
catchments, since the selected GCMs provide both an increase and a 
decrease in future precipitation. It is also noted that the uncertainty 
associated with the NCV is an important source of uncertainty depending 


on the emission scenario considered. 


It is well described in the literature that the largest source of 
uncertainty in climate change impact studies on water resources is the 
choice of the GCM (e.g., Graham, Hagemann, Jaun, € Beniston, 2007). 
Accordingly, our results show that climate model uncertainty is a key 
factor for impact studies. However, the uncertainty associated with the 
NCV cannot be neglected in hydrological climate change impact studies. 
In the present work, the results demonstrate that the uncertainty 
related to the NCV is critical in the evaluation of the impacts of climate 
change on high flows, and to some extent, comparable to the GCM 


uncertainty and emission scenario uncertainty. 


Future work should consider additional basins in order to explore 
the contrasted climate types and climate variability that determine the 
hydrological regimes of the Mexican territory. In addition, other sources 
of uncertainty should be evaluated, such as the bias correction method 
and the hydrological model uncertainties. Finally, the present work 
based on the CMIP3 multi-model dataset will be next validated 
according to the recent climate model simulations from the CMIP5 multi- 


model dataset. 
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